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Motivated by recent Raman and resonant inelastic X-ray 
scattering experiments performed for Mott insulators, which 
suggest formation of excitons in these systems, we present a 
theory of exciton formation in the upper Hubbard band. The 
analysis based on the spin polaron approach is performed in 
the framework of an effective t-J model for the subspace of 
states with one doubly occupied site. Our results confirm the 
existence of excitons and bear qualitative resemblance to ex- 
perimental data despite some simplifications in our approach. 
They prove that the basic underlying mechanism of exciton 
formation is the same as that which gives rise to binding of 
holes in weakly doped antiferromagnets. 

PACS numbers: 74.25.Jb, 71.10.Fd, 71.27.-Fa 



The electronic properties of Mott insulators, like high 
temperature superconductors in the limit of low doping, 
belong to the most intriguing questions which faces con- 
temporary condensed matter physics. Properties of such 
systems are determined by the large onsite Coulomb in- 
teraction U or the charge transfer energy A of removing 
an electron from the p orbital of oxygen and putting it 
into the d orbital of the transition metal. The struc- 
ture of occupied electronic states has been intensively in- 
vestigated by means of the angle-resolved photoemission 
spectroscopy (ARPES) and much is already known about 
the electronic structure below the insulating gap. Natu- 
ral probes of unoccupied states above the gap are Raman 
[01 and resonant inelastic X-ray scattering Q (RIXS). 

The Raman scattering provides information about ex- 
citations at k = 0. If we ignore polaritonic effects the 
cross section for the Raman scattering may be obtained 
directly from the golden rule. The final formula con- 
tains matrix elements evaluated between an ini- 
tial state i and a final state /. In the case of the Mott 
insulator in two dimensions (2D) described by the Hub- 
bard model at half-filling in the limit of strong correla- 
tions (large U), which is a natural framework to discuss 
such a system, the initial state i is just the ground state of 
the quantum antiferromagnet (AF) in 2D, while the final 
state may be identified as an excitation in the AF state. 
The electronic states \n) above the insulating gap appear 
at the intermediate stage of the Raman process, which 
may be seen from the explicit formula for {f\Mfi\i), 
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Ci and e/ are energies of the initial and final states, while 
huji, Bi and hiVf, ej represent energies and polarizations 
of incident and outgoing photons. The current operator 
for the Hubbard model on the square lattice is given by, 
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The summation is carried over lattice sites i, vectors 5 
which connect each site with its four nearest neighbors 
and spin polarizations a. The Hubbard model in the 
lower Hubbard band at half-filling in the limit of strong 
correlations is equivalent to the Heisenberg model which 
action is restricted to the Hilbert space spanning states 
that correspond to configurations in which each site is 
occupied by a single spin-1/2 fermion. The current op- 
erator (|^) acting on the antiferromagnetic state, shifts a 
fermion to a nearest neighbor site that is already occu- 
pied, which means that it couples lower and upper Hub- 
bard bands (LHB, UHB). The newly created hole and 
the double occupancy (DO) of a site reside on nearest 
neighbor sites. We conclude, that in order to contribute 
some substantial weight to the Raman spectra, the in- 
termediate state should be a bound state of some sort, 
which is formed by the hole and the DO that plays the 
role of a particle. In other words, we expect formation 
of excitons in the UHB. The current operator at k = 
is a vector and obeys a dipole selection rule. Since the 
ground state of the quantum AF is fully symmetric, the 
current operator couples it only with excitons of p-wave 
symmetry. 

The RIXS process may be interpreted in a slightly sim- 
ilar way, which in addition requires explicit consideration 
of the orbital structure of involved states. According to 
a most likely scenario for the undoped cuprates 
3d electrons interact with the Is core hole created during 
the dipole transition of an Is electron to the 4p orbital 
which accompanies the absorption of an incident photon. 
This strong interaction gives rise to a shift of the 3d elec- 
trons across the gap. In the last step, the Ap electron 
returns to the Is orbital and an outgoing photon is emit- 
ted. The RIXS process like the Raman scattering creates 
a particle-hole pair in the UHB which carries energy and 
in addition has non-vanishing momentum. It is easy to 
imagine three dispersing modes in the UHB with a sin- 
gle hole which might be relevant to the RIXS physics. 
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a positively charged spin polaron related with the hole, 
a negatively charged spin polaron related with the DO 
and a bound state of these two entities (exciton). Our 
calculation will clarify the mechanism of formation and 
propagation of such objects. 

Since we are concerned with the the Hubbard model 
in the large-J7 limit we choose a standard way of sepa- 
rating Hubbard bands by means of a unitary transfor- 
mation. A version of an effective t-J model obtained in 
this way, which acts in the UHB was derived by Eskes 
et al. Q . Physical processes contained in that model are 
actually very similar to the processes which are present 
in the standard t-J model defined in the LHB, but in- 
volvement of doubly occupied sites prohibits derivation 
of a compact form for the Hamiltonian in the case of the 
UHB. The effective model for the UHB promotes AF cor- 
relations, like its counterpart for the LHB, and therefore 
we choose the Neel state |iV) as a point of reference by 
means of which we are going to present the Hamiltonian 
in a comprehensive way. The virtue of this representa- 
tion is that it demonstrates the crucial duality between 
a hole and a DO. For the purpose of this paper we tradi- 
tionally neglect in the effective Hamiltonian terms related 
with correlated hopping if they do not involve simultane- 
ously a hole and a DO. Let a{i) denote the direction of 
spin at the site i in the Neel state |A^) and d-{i) = —cr{i). 
States spanned by the fermionic Hilbert space may be ob- 
tained by acting on |7V) with products of ct|.^^Ccr(i) = 

^ d\. They invert a spin (create a 



-cr(i) 



magnon), create a hole or a DO at the site i. The rep- 
resentation in terms of products will be unique if each 
index z, which refers a to site, appears in all products 
no more than once. Since the effective Hamiltonian ob- 
tained by Eskes et al. corresponds to the second-order 
perturbation expansion it acts maximally on three sites 
which lie in row in the lattice and hence affects only three 
operators in the above mentioned products. For the same 
reason only a single hole-creation operator and a single 
DO-creation operator may be involved in the action of 
the effective Hamiltonian. In (^^5|) the action of the t-J 
model in the UHB with a single hole is represented by 
means of states which are coupled by this operator. Re- 
lations d^p^ actually demonstrate how the action of the 
effective Hamiltonian interchanges pairs or trios of cre- 
ation operators represented by O and O' between prod- 
uct representations like AOB\N) and AO' B\N) of these 
states, [A and B denote products of operators which act 
on different sites than those involved in the action of O 
and O'). Lack of an operator for a given site i in the 
product is formally represented by the identity operator 
li. In this case the site i is occupied by a single spin 
which points in the direction determined by the original 
Neel configuration. The value of the matrix element of 
the Hamiltonian is given in (^[]l5|) after a colon. The 
energy of the Neel state is chosen as a point of refer- 



ence from which we count the diagonal contributions to 
the Hamiltonian. A different description of relations 
[Tsl ) is that they represent a sum over states with factors 
given by amplitudes presented after colons. That sum 
is created by action of the Hamiltonian on a given state 
represented by a product. Each component of the sum 
may be obtained by exchanging operators on two or three 
sites according to the rules represented by (|-15). Pro- 
cesses in (||) represent hopping, or in other words shifts 
of a hole or a DO accompanied by creation or annihi- 
lation of a magnon. In anti-parallel spins on a pair 
of nearest sites are simultaneously inverted. (^[]^) are 
related with the rise in the contribution to diagonal ma- 
trix elements of the Hamiltonian (which play the role of 
potential energy), in comparison to the vanishing contri- 
butions from (^), which appears if nearest sites are not 
occupied by parallel spins. That rise is higher if such a 
pair is occupied by a hole and a DO, which may be seen 
in the second term in (||). On the other hand in such 
a case an additional process is possible that may lower 
the energy, namely the exchange of the hole and the DO 
which may be seen in (^). Processes described by formu- 
las ([lO|-p^ represent exchange, rotations and shifts of the 
hole-DO pair between nearest links which in some cases 
are accompanied by a shift of a magnon. 
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The duality between the hole and the DO may be clearly 
seen in (||-|l5|). We may also notice the similarity with the 
t-J model acting in LHB. The main differences are that 
the rise of the potential energy is only J/2 if a hole pair 
occupies nearest sites and not J like for a holc-DO pair 
and that in the LHB only one kind of charged particles 
(holes) is present. Nevertheless, one may expect that the 
possible mechanism of binding of a hole and a DO in the 
UHB is similar to the mechanism of hole pairing in the 
LHB. 

Much is already known about binding of holes in 
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weakly doped antiferromagnets |10|. We will use these 
insights to analyze formation of an exciton in the UHB. 
Fast motion, with the rate ^ t, mediated by terms (|^) 
of a hole or a DO created at some site in the Neel state 
produces magnons which lie on the track of the hole or 
the DO and contribute to the rise in the potential energy 
(terms (^[^)). Strings consisting of magnons are pinned 
to the site where these objects have been initially created 
and form a potential well for a hole and a DO which con- 
fines their motion. Much slower processes related to the 
process of inverting a pair of anti-parallel spins at near- 
est neighbor sites (process (^)) at the begin of the string, 
which occur at the rate ~ J shorten strings and give rise 
to coherent propagation of both objects. Since the sep- 
aration of energy scales for hopping of a hole or a DO 
and annihilation of defects (magnons) is pronounced it is 
plausible to introduce the notion of a spin polaron. We 
may define in this context a hole or a DO spin-polaron 
as a solution of the Schrodinger equation for a particle 
in the potential well, the wave function of which is as a 
combination of states that may be reached by hopping of 
a hole (or a DO) created at an initial site i, 
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Vi parameterizes the geometry of a path along which the 
hole or a DO has been moving and \V^''^'~'^) is a state 
which has been created in this way. To be precise, at this 
stage of our considerations we analyze an " unperturbed" 
Hamiltonian which consists of terms symbolized by re- 
lations (^) and (H^) and represents fast oscillations of a 
hole or a DO in the vicinity of the initial site. Slower pro- 
cesses are neglected in this phase of analysis. "Orbital" 
states of polarons created at all possible sites exhaust, 
in principle, the relevant portion of the Hilbert space. In 
practice, the calculation may be confined to some low ex- 
cited states or even to the polaron ground state. The con- 
cept of spin polarons may be extended to the description 
of interaction between a hole and a DO. If we consider 
the hole and the DO which have been initially created at 
distant sites we just apply the generalization of a single 
particle problem by assuming that the wave function in 
this case is a product of wave functions for a hole and a 
DO. For a pair of particles created an nearest-neighbor 
sites we define a hole-DO spin-bipolaron I'^fi^f) as 
a sum over states which may be reached by independent 
hopping of a hole and a DO created at the sites i and j, 
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The meaning of symbols is the same as in the case of 
single particles. At this stage of considerations we pro- 
hibit, by definition, the hole and the DO to follow along 
the trace left by the accompanying particle and solve the 
Schrodinger equation in this restricted Hilbert space for 



an "unperturbed" Hamiltonian, which takes into account 
only hopping of the hole or the DO and the diagonal " po- 
tential" contribution to the energy from the magnetic J- 
term in the t-J model. By means of that restriction we 
achieve that spin bipolarons are localized and we may 
proceed as in the case of a single hole. We make a fur- 
ther approximation and neglect some path details includ- 
ing the possibility of path crossing. It is self-evident that 
the coefficients a-p- and a-p.^-p. for polarons in the ground 
state will depend only on the lengths /j,, v of paths Vi and 
Vj, and thus a-p^ = and aVi,Vj = ot^.y By calculat- 
ing the matrix elements of the full Hamiltonian in the po- 
laron basis we derive an effective Hamiltonian expressed 
in the polaron language. Some processes, like oscillations 
of the hole or the DO in the vicinity of a polaron center, 
are already incorporated into the eigenenergy of polaron 
states. The rest which we discuss now either renormalizes 
the eigenenergies of polarons or gives rise to off-diagonal 
matrix elements in the polaron Hamiltonian. The latter 
circumstance occurs in the case of shortening of strings 
attached to a single particle. This process governs propa- 
gation of a single hole or a single DO in an antiferromag- 
netic spin background. Formation of bipolarons, and, in 
particular, strings that connect a hole and a DO, is an 
effective way of lowering the energy. A compromise be- 
tween two opposite tendencies to minimize the kinetic 
energy of holes and DOs and to reduce the disturbance 
of the antiferromagnetic background is reached in this 
way. By shrinking at one end and expanding at the op- 
posite end, a string may move, while keeping a moderate 
length. The physics of interaction between a hole and 
a DO is richer than just for two holes, because it incor- 
porates processes described by relations (|9|-p^, but the 
method of deriving the Hamiltonian in the pol aron lan- 
guage is the same as for one kind of particles . Instead 
of dwelling on details we present now an explicit form of 
the effective Hamiltonian H^ff- h\ and d\ are fermionic 
operators which create a hole and the DO spin-polaron 
at the site i. They obey the relation h\d\ — Q because 
none of sites may be simultaneously occupied by a hole 
and a DO. The bipolaron j)*^) is created by the pair 

where i and j are nearest- neighbor sites. 

Heff^ El Y.^h% + d\di) + 2X1 ^1+5+5'^^ 

■i i,6,S' -.S'^ — S 

+ 2x1 Y d\^,^,,d, + {E2-2Ei) 

i,S,S':S'^~5 

X Y 4+5di+sh\hi + {E2UJ1 + ri + i/2) 

1,5 

X Y {4di+5'h\^ghi + H.c.) + {E2UJ2 + T2) 

i.5,S':5'^5 

X E idl+s+5'd'M+s+5'+S"h^+S + H.C.) 

i,SJ',S":S'yi-S,S"yi-S' 
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i,S i,S.5':S'^6 

- 2x1 dl^gdi+shl^g^g,hi 

+ X2 {d\+s+snd,+sh\^g,h, + H.c.). 

(18) 

El and denote the eigenenergies of the spin polaron 
and the bipolaron which represent the energy of a particle 
or two particles pinned to starting-point sites by strings, 
while the rest of parameters is either related with overlap 
between polaron states at different sites or processes that 
couple polaron states: 

^i^ Y (2-l)''+""'a,,^a.-i.p+i, (19) 

Y (^-l)''+""'a.,Ma>'-2,M+2, (20) 

n = t ^(z - l)^-ia^,oa^-i,o, (21) 

ra =t^(z- l)^-2^^^^^^_^^^^ (22) 

Xi = J/2^(z-l)^-2^^^^_^^ (23) 

X2 = J/2 Y (^-l)^+""'a.,.aM-i«'^-i- (24) 

t = J/2(ao,o)'. (25) 

As we have mentioned before polaron states are not nec- 
essarily orthogonal and in order to find eigenenergies of 
the system one has to diagonalize the operator 0~j^jHf,f f , 
where the operator Oe/ / represents overlap between po- 
laron states, 

Oeff= l+uji Y, {dl^i+S'hl^s^i + H.c.) +UJ2 
^ Y i4+s+s'd-A+s+s'+s"^^+s + H.c.). 

i,S,S',S":S'^-S,S"^-S' 

(26) 

The Hamiltonian (|l^) and the overlap operator (^6|) for 
a hole and a DO are very similar to operators that repre- 
sent interaction of holes in the LHB. Thus the underlying 
mechanisms of exciton formation, if it occurs, and pairing 
of holes should be very similar. Fig.|l| depicts the energy 
dispersion of eight states with lowest energies obtained 
by solving the eigenvalue problem for the operators (llq) 



and (^6|) in the case of a single hole and a single DO. The 
line in the upper part of Fig.|l| represents the lower bound- 
ary for the energy of scattering states. There exist three 
bound excitonic states at the wave vector k = (0, 0). The 
ground state has ti-wave symmetry. The pair of bound 
states with higher energy has p-wave symmetry, which 
guarantees that they may take part in dipole transitions 
and contribute to Raman spectra. 
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FIG. 1. Energy dispersion of eight lowest states in the UHB 
with a single hole on a cluster of 32 x 32 sites for J/t = 0.4. 

For other values of the wave vector up to five exciton 
states at different energies may exist. The typical dis- 
tance between the hole spin-polaron and the DO spin- 
polaron in the exciton state is one lattice spacing which 
indicates that they are tightly bound and constitute a 
bipolaron. Among dispersion curves which form in Fig.|l| 
one may identify bands which shape agree with the form 
of energy-dispersion curves for the particle-hole excita- 
tion measured by Hasan et al. in the RIXS experiment, 
but the bandwidth in the theoretical approach for a rea- 
sonable hopping parameter t = 0.35 eV seems to be much 
to small, which should be attributed to the fact that we 
have neglected the possibility of direct hopping to further 
neighbors. Additional analysis is also needed to under- 
stand the dependence of spectral intensity on the position 
of the wave vector in the Brillouin zone. 
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